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Abstract 

We have determined the abundance of i'-band drop-outs in the HST/ ACS GOODS 
surveys and the Hubble Ultra Deep Field (UDF). The majority of these sources 
are likely to be z ~ 6 galaxies whose flux decrement between the F775W i'-band 
and F850LP z'-band arises from Lyman-alpha absorption. We have shown with 
Keck/DEIMOS and Gemini/GMOS spectroscopy that this technique does indeed 
select high redshift galaxies, and we discovered Lyman-a emission in the expected 
redshift range for about a third of the galaxies with z' AB < 25.6 in the 150arcmin 2 
of the GOODS-South field. The i-drop number counts in the GOODS-North field 
are consistent, so cosmic variance is possibly not the dominant uncertainty. The 
increased depth of UDF enables us to reach a ~ 10 a limiting magnitude of z' AB = 
28.5 (equivalent to 1.5 hj^ Mq yr _1 at z = 6.1, or 0.1 L^y for the z ~ 3 [/-drop 
population). The star formation rate at z ~ 6 was approximately x6 less than 
at z ~ 3. This declining comoving star formation rate (0.005 hjQ Mq yr _1 Mpc -3 
at z ~ 6 at Luy > 0.1 L* for a Salpeter IMF) poses an interesting challenge for 
models which suggest that Ljjy > 0.1 L* star forming galaxies at z ~ 6 reionized the 
universe. The short-fall in ionizing photons might be alleviated by galaxies fainter 
than our limit, or a radically different IMF. Alternatively, the bulk of reionization 
might have occurred at z ^> 6. We have recently discovered evidence of an early 
epoch of star formation in some of the z'-drops at z ~ 6. Spitzer images with IRAC at 
3.6 — 4.5 /um show evidence of the age-sensitive Balmer/4000 A, dominated by stars 
older than 100 Myr (and most probably 400 Myr old). This pushes the formation 
epoch for these galaxies to Zf OTm = 7.5 — 13.5. There are at least some galaxies 
already assembled with stellar masses ~ 3 x 10 10 M Q (equivalent to 0.2 M* today) 
within the first billion years. The early formation of such systems may have played 
a key role in reionizing the Universe at z ~ 10. 
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1 Introduction 



There has been enormous progress over the past decade in discovering galaxies 
and QSOs at increasingly high redshifts. We are now probing far enough back 
in time that the Universe at these early epochs was fundamentally different 
from its predominantly ionized state today. Observations of z > 6.2 QSOs 
(Becker et al. 2001, Fan et al. 2002) show near-complete absorption of flux 
at wavelengths short- ward of Lyman-a (Gunn & Peterson 1965), indicating 
that the Universe is optically thick to this line, and that the neutral fraction 
of hydrogen is much greater than at lower redshifts. Initial cosmic microwave 
background results from the WMAP satellite indicated that the Universe was 
completely neutral at redshifts of z ~ 10 (Kogut et al. 2003). There is an 
ongoing debate as to what reionized the Universe at z > 6: is it AGN or ion- 
izing photons (produced in hot, massive, short-lived stars) escaping from star 
forming galaxies? Along with the escape fraction for these ionizing photons, 
it is crucial to know the global star formation rate at this epoch. 

Here we detail our work in identifying star-forming galaxies within the first 
billion years of the Big Bang through deep imaging with the Hubble Space tele- 
scope (HST). Our selection method is based on the Lyman break technique, 
and we have shown this is effective in isolating z ~ 6 by obtaining spectro- 
scopic confirmation with the 10-m Keck telescopes. Our analysis of the Hubble 
Ultra Deep Field (the deepest images ever obtained) enables us to explore 
the faint end of the luminosity function (the bulk of galaxies with low star 
formation rates). Our discovery of this i'-drop galaxy population is used to 
infer the global star formation rate density at this epoch (z « 6). We consider 
the implications for reionization of the UV flux from these galaxies. We also 
use infrared data from the Spitzer Space Telescope to determine the spectral 
energy distributions (SEDs) of our i'-drops, and constrain the previous star 
formation histories of these sources. 

Throughout we adopt the standard "concordance" cosmology of Vt M = 0.3, 
VIa = 0.7, and use h 70 = iJ /70kms _1 Mpc -1 . All magnitudes are on the AB 
system. 



2 The Lyman-Break Technique at z ~ 6 

In order to select z ~ 6 galaxies, we use the "Lyman break technique" pi- 
oneered at z ~ 3 using ground-based telescopes by Steidel and co-workers 
(Steidel, Pettini & Hamilton 1995; Steidel et al. 1996) and using HST by 
Madau et al. (1996). At z ~ 6, we can efficiently use only two filters, above 
(z'-band 9000A) and below (i'-band 8000A) the continuum break at Lyman-a 
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Fig. 1. Left: The ACS-z' and -z' bandpasses overplotted on the spectrum of a 
generic z = 6 galaxy (solid line), illustrating the utility of our two-filter technique 
for locating z ~ 6 sources. Right: Model colour-redshift tracks for galaxies with 
non-evolving stellar populations (from the template spectra of Coleman, Wu & 
Weedman 1980, ApJS 43, 393 ). The contaminating 'hump' in the (i' — z') colour at 
z ps 1 — 2 arises when the Balmer break and/or the 4000 A break redshifts beyond 
the i'-filter. 

line (A = (1 + z) x 1216 A) since the integrated optical depth of the Lyman-o; 
forest is ^ 1 (see Figure 1) and there is essentially no flux at shorter wave- 
lengths. The key issue is to work at a sufficiently-high signal-to-noise ratio that 
Lyman break "drop-out" galaxies can be safely identified through detection 
in a single redder band (i.e., z'-band). This approach has been demonstrated 
to be effective by the SDSS collaboration in the detection of z ~ 6 QSOs 
using the i'- and z'-bands alone (Fan et al. 2002), using wide- area but shal- 
low ground-based imaging (sensitive to bright but rare quasars). Our group 
has analysed deep imaging from the Advanced Camera for Surveys (ACS) 
on HST with the same sharp-sided SDSS F775W (i') and F850LP (z') niters 
to locate the"i-drops" which are candidate z ~ 6 galaxies (fainter but more 
numerous than the QSOs). In Figures 1&2 we illustrate how a colour cut of 
(i' — z')ab > 1-3 can be effective in selecting sources with z > 5.6. 

We have analysed the ACS images from the HST Treasury "Great Observatory 
Origins Deep Survey" (GOODS; Giavalisco & Dickinson 2003) to discover this 
i'-drop population of z ~ 6 galaxies. In Stanway, Bunker & McMahon (2003) 
we searched 150arcmin 2 of the GOODS-South field (the Chandra Deep Field- 
South) to select (i' — z') AB > 1.5 objects, of which six were probable z > 5.7 
galaxies brighter than zab < 25.6. Subsequent papers by the GOODS team 
(Dickinson et al. 2004; Giavalisco et al. 2004) reproduced many of these %'- 
drop candidates. To address potential cosmic variance issues, we performed a 
similar analysis in the GOODS-North field, which yielded a consistent estimate 
of the surface density of z ~ 6 star forming sources (Stanway et al. 2004b). 
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2.1 Spectroscopic Confirmation of z ps 6 Galaxies 



The effectiveness of the Lyman-break technique has been demonstrated at 
lower redshifts through spectroscopic confirmation of hundreds of [/-band 
dropouts at z ~ 3 and £?-band drop-outs at z ~ 4. However, it is impor- 
tant to demonstrate that the i'-drop selection is similarly isolating galaxies at 
z ~ 6, if we are to use the surface density of ? '-drops to draw global inferences 
about galaxies within the first billion years. A known possible contaminant is 
the Extremely Red Object (ERO) population of evolved galaxies at z ~ 1 — 2 
which can produce large {%' — z') colours: deep near-infrared imaging (e.g., 
from the i/ST/NICMOS survey of the UDF) should identify EROs, and we 
consider this possible contamination Stanway, McMahon & Bunker (2005). 
Low- mass Galactic stars (M/L/T dwarfs) are another interloper population. 



We have obtained deep spectroscopy with Keck/DEIMOS of several of the 
brightest i'-drops identified by Stanway, Bunker & McMahon (2003) in GOODS- 
South; in the case of SBM03#3 {z! = 25.6) and , an i'-drop in the Chandra 
Deep Field South photometrically-selected from HST/ACS images to lie at 
z 6. 



For two of these objects, we discovered single emission lines (Fig. 3) at 8245 A 
& 8305 A, consistent with Lyman-a emission from galaxies at z — 5.78 (SBM03#3, 
Bunker et al. 2003) and z = 5.83 (SBM03#1, Stanway et al. 2004a). The 
spectrally-resolved profiles of the emission lines are asymmetric (as high-z 
Lyman-a tends to be) with a P-Cygni-like profile and a sharp cut-off on 
the blue wing (Fig. 2). The line fluxes are ps 2 x 10 _17 ergscm _1 s _1 . The 
equivalent widths are W/ r e S t = 20 - 30 A using the z'-band photometry from 
HST/ACS, which is within the range seen in high- z star-forming galaxies 
at z ~ 3 — 4. The velocity widths of the Lyman-a are comparatively nar- 
row (vfwhm = 300 km s" 1 ) and we do not detect the high-ionization line 
N v A 12140 A, which strongly support the view that this line emission is pow- 
ered by star formation rather than an AGN. Our spectroscopic redshifts for 
these objects confirms the validity of the i'-drop selection technique to select 
star-forming galaxies at z ~ 6. Subsequently, a programme of spectroscopy 
with GMOS "nod & shuffle" on Gemini (GLARE, Stanway et al. 2004b) has 
confirmed more galaxies in GOODS-South. Slitless spectroscopy of the UDF 
with HST/ACS (GRAPES, Malhotra et al. 2005) has also confirmed the spec- 
tral breaks of the i'-drop galaxies, although the spectral resolution of R ~ 70 
is insufficient to see the Lyman-a emission lines. 
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Fig. 2. Left: Colour-magnitude diagram for the UDF data with the limit z' AB < 28.5 
and (i' — z')ab = 1-3 colour cut shown (dashed lines). Such a catalogue could be 
contaminated by cool stars, EROs and wrongly identified extended objects and 
diffraction spikes but nonetheless provides a secure upper limit to the abundance 
of z« 6 star forming galaxies. Circles and arrows (lower limits) indicate our i'-drop 
candidate z « 6 galaxies. The asterisk is the only unresolved i'-drop in our UDF 
sample, a probable star. The solid circle is the brightest i'-drop in the UDF, the 
spectroscopically-confirmed galaxy SBM03#1 at z = 5.83, with our discovery spec- 
trum from Keck/DEIMOS shown (right, Stanway et al. 2004b. 
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Fig. 3. Our 2D Keck/DEIMOS spectrum confirming the brightest i'-drop in the 
UDF is indeed az«6 galaxy: prominent Lyman-a emission is seen at z = 5.83 
(from Stanway et al. 2004b). There has been independent spectroscopic confirmation 
by Dickinson et al. (2004). 
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Fig. 4. Left: Cumulative source counts per arcmin 2 of i'-dropouts as a func- 
tion of z'-band magnitude. The new UDF data (over an area of 11 arcmin 2 for 
28.5 > z' AB > 27.0) is compared with our z' AB < 25.6 single epoch GOODSv0.5 
ACS/WFC imaging over 300 arcmin 2 (Stanway, Bunker & McMahon 2003; Stan- 
way et al. 2004), and combined 5 epoch GOODSvl.O images to z' AB < 27.0 (Stanway 
2004, PhD. thesis, Cambridge University). Also plotted are the surface densities of 
i'-drops from Giavalisco et al. (2004) Bouwens et al. (2003, 2004) 
Right: An updated version of the 'Madau-Lilly' diagram (Madau et al. 1996; Lilly 
et al. 1996) illustrating the evolution of the comoving volume- averaged star forma- 
tion rate. Our work from the UDF data is plotted a solid symbol. Other determi- 
nations have been recalculated for our cosmology and limiting UV luminosity of 
1.5 hjQ M & yr -1 at z = 6.1 (equivalent to 0.1 L3 at z 3 from Steidel et al. 1999), 
assuming a slope of a = —1.6 for z > 2 and a = —1.3 for z < 2. Data from the 
CFRS survey of Lilly et al. (1996) are shown as open circles; data from Connolly et 
al. (1997) are squares; and the Lyman break galaxy work of Steidel et al. (1999) is 
plotted as crosses, of Fontana et al. (2002) as inverted triangles and that by Iwata et 
al. (2003) as an open diamond. Pentagons are from Bouwens, Broadhurst & Illing- 
worth (2003), upright triangles are the GOODS i'-drop results from Giavalisco et 
al. (2004). The three ACS estimates of Bouwens et al. (2003) are shown by small 
crossed circles and indicate three different completeness corrections for one sample 
of objects - the larger symbol is the recent re-determination using a new catalogue 
by this group from a deeper dataset (the UDF flanking fields - Bouwens et al. 2004). 



3 The star formation rate density at z ~ 6 



We base our measurement of the star formation rate on the intensity of the 
rest-frame UV continuum emission, redshifted into the z'-band at z ~ 6; this 
is dominated by the light from the most massive stars (M > 10 M Q which are 
the hottest and bluest stars), which are the shortest lived, and hence a tracer 
of the instantaneous star formation rate. To go from the number of massive 
stars formed to the total star formation rate requires an assumption about 
the stellar initial mass function (IMF). In the absence of dust obscuration, the 
relation between the flux density in the rest-UV around 1500 A and the star 
formation rate (SFR in M Q yr -1 ) is given by L uv = 8 x 10 27 SFRergss _1 Hz -1 



6 



from Madau, Pozzetti & Dickinson (1998) for a Salpeter (1955) IMF with 
0.1 M < M* < 125 M . Our relatively bright magnitude cut of z' AB < 25.6 
for the GOODSv0.5 data (Stanway, Bunker & McMahon 2003; Stanway et al. 
2004a) corresponds to an unobscured star formation rate of 15 /if 2 M Q yr" 1 at 
z = 5.9 (the luminosity- weighted average redshift), equivalent to L* uv for the 
[/-band dropout population at z ~ 3. 

We recover the z ~ 6 rest-frame UV luminosity function from the observed 
number counts of i'-drops to faint magnitudes in the UDF. Although our colour 
cut selects galaxies with redshifts in the range 5.6 < z < 7.0, an increasing 
fraction of the z'-band flux is attenuated by the redshifted Lyman-a forest. 
At higher redshifts we probe increasingly shortward of A rcs t = 1500 A (where 
the luminosity function is calculated) so the /^-corrections become significant 
beyond z ~ 6.5. 

We now compare the comoving star formation rate deduced for z ~ 6 galaxies 
based on our candidate i'-dropout source counts with predictions based on a 
range of rest-frame UV luminosity functions. For convenience we assume that 
there is no evolution over the sampled redshift range, 5.6 < z < 6.5, spanned 
by our i'-drop selection (equivalent to a range between 0.8 — 1.0 hj^ Gyr after 
the Big Bang). We take as a starting point the luminosity function for the well- 
studied Lyman-break [/-dropout population, reported in Steidel et al. (1999), 
which has a characteristic rest-UV luminosity Lg = 15 /if 2 M yr _1 and a char- 
acteristic comoving number density at z ~ 3 is = 0.0014 /i^Mpc -3 mag" 1 
in our cosmology. The faint end slope of the Schechter function at z ~ 3 is 
relatively steep (a = —1.6) compared with a = —1.0 to —1.3 for lower-redshift 
galaxy samples (e.g., Blanton et al. 2003). 

A tantalizing result from our work so far is that at z ~ 6 there are far fewer 
UV-luminous star forming galaxies than would have been predicted if there 
was no evolution, based on a comparison to the well-studied z ~ 3 — 4 Lyman 
break population (Steidel et al. 1999). In fact, the volume averaged (comoving) 
star formation rate in galaxies with > 15 hjQ M Q yr" 1 is ~ 6x less at z ~ 6 
than at z ~ 3. So it would appear that at this crucial epoch, where abundant 
sources of UV photons are needed to reionize the Universe, there is in fact 
a deficit of luminous star forming galaxies. We have conclusively shown that 
the bright end of luminosity function has evolved greatly from z ~ 6 to z ~ 
3. Other groups have claimed less dramatic evolution or even no evolution 
in the total volume-averaged star formation rate, based on the same fields 
(Giavalisco et al. 2004; Dickinson et al. 2004) and similar HST/ ACS data sets 
(Bouwens et al. 2003; Yan, Windhorst & Cohen 2003), but these groups work 
closer to the detection limit of the images and introduce large completeness 
corrections for the faint source counts. To do a complete inventory of the UV 
light from star formation we must address the contribution of low-luminosity 
star-forming galaxies to the ionizing flux. The public availability of the Hubble 
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Ultra Deep Field (UDF; Beckwith, Somerville & Stiavelli) can address this 
puzzle by pushing down the luminosity function at z ~ 6 to well below the 
equivalent of L* for the z ~ 3 population. 



4 The Hubble Ultra Deep Field 



The Hubble Ultra Deep Field (UDF) was a Cycle 12 STScI Director's Dis- 
cretionary Time programme executed over September 2003 - January 200, 
comprising 400 orbits in 4 broad-band filters (including F775W i'-band for 
144 orbits; F850LP z'-band for 144 orbits). The UDF lies within the Chandra 
Deep Field South (CDF-S). As the UDF represents the deepest set of images 
yet taken, significantly deeper than the /-band exposures of the Hubble Deep 
Fields (Williams et al. 1996), and adds the longer- wavelength z'-band, it is 
uniquely suited to the goals of our program. 

We performed the first analysis of i'-drops in the UDF (Bunker et al. 2004), 
presenting details of 54 candidate star forming galaxies at z ~ 6 in a preprint 
the day after the release of the UDF data. Our analysis was subsequently 
independently repeated by Yan & Windhorst (2004), with agreement at the 
98% level (see Bunker & Stanway 2004 for a comparison). We take our mag- 
nitude limit as z' AB < 28.5 (a 10 o cut). We measure a star formation den- 
sity of 0.005 h-jQ Mq yr _1 Mpc -3 at z ~ 6 from galaxies in the UDF with 
SFRs > 1.5 hjQ M yr _1 (equivalent to 0.1 L* uv at z « 3). 

At the relatively bright cut of z' AB < 25.6 used in Stanway, Bunker & McMa- 
hon (2003) from the GOODSv0.5 survey, the UDF data is 98% complete for 
sources as extended as = 0.5arcsec. Interestingly, we detect no extended 
(low surface brightness) i'-drops to this magnitude limit in addition to the tar- 
getted i'-drop SBM03#1 in the deeper UDF data. This supports our assertion 
that the i'-drop population is predominantly compact and there cannot be a 
large completeness correction arising from extended objects (c.f. Lanzetta et 
al. 2002). The ACS imaging is of course picking out Hn star forming regions, 
and these UV-bright knots of star formation are typically compact (< lkpc, 
< 0.2arcsec at z « 6) even within large galaxies at low redshift. 

It is interesting that the level of stellar contamination in the UDF i'-drops is 
only 2%, compared with about one in three at the bright end (z' AB < 25.6; 
Stanway, Bunker & McMahon 2003; Stanway et al. 2004). This may be because 
we are seeing through the Galactic disk at these faint limiting magnitudes to 
a regime where there are no stars (see also Pirzkal et al. 2005). 
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5 Implications for Reionization 



We compare our i'-drop luminosity function with the work of Madau, Haardt & 
Rees (1999) for the density of star formation. We have updated their equation 
27 for the more recent concordance cosmology estimate of the baryon density 
of Spergel et al. (2003), and for the predicted mean redshift of our sample 
(z = 6.0): 



This relation is based on the same Salpeter Initial Mass Function as we have 
used in deriving our volume-averaged star formation rate. C is the concen- 
tration factor of neutral hydrogen, C = (phi) (phi)~ 2 - Simulations suggest 
C ~ 30 (Gnedin & Ostriker 1997). The escape fraction of ionizing photons 
(/esc) fc> r high- redshift galaxies is highly uncertain (e.g., Steidel, Pettini & 
Adelberger 2001), but even if we take f esc = 1 (no absorption by H i) this esti- 
mate of the star formation density required is a factor of m 5 higher than our 
measured star formation density of 0.005 h 70 M Q yr~ x Mpc^ 3 at z ~ 6 from 
galaxies in the UDF with SFRs> 1.5 /if 2 M yr -1 . For faint end slopes of 
a — 1.8 — > —1.3 galaxies with L > 0.1 L* account for 32 — 80% of the total 
luminosity, so would fall short of the required density of Lyman continuum 
photons required to reionize the Universe. If the faint-end slope is as steep as 
a ~ —1.9 then there would just be enough UV Lyman continuum photons 
generated in star forming galaxies at z ~ 6 (assuming a Salpeter IMF), but 
the required escape fraction for complete reionization would still have to be 
implausibly high (f esc ~ 1, whereas all high- ,2 measurements to date indicate 
f esc <C 0.5: Fernanadez-Soto, Lanzetta & Chen 2003; Steidel, Adelberger & 
Pettini 2001). 

AGN are also under-abundant at these epochs (e.g., Dijkestra, Haiman & Loeb 
2004). If star forming galaxies at redshifts close to z = 6 were responsible for 
the bulk of reionization, then a very different initial mass function would be 
required, or the calculations of the clumping factor of neutral gas would have 
to be significantly over-estimated (see also Stiavelli, Fall & Panagia 2004). Al- 
ternatively another low- luminosity population (e.g., forming globular clusters; 
Ricotti 2002) could be invoked to provide some of the shortfall in ionizing 
photons. It is also plausible that the bulk of reionization occurred at redshifts 
well beyond z = 6: the WMAP polarization data indicate z re i on > 10 (Kogut 
et al. 2003), and it is possible that the Gunn- Peterson troughs seen in z > 6.2 
QSOs (Becker et al. 2001; Fan et al. 2002) mark the very last period of a 
neutral IGM. 
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Fig. 5. Left: The Spitzer IRAC and HST ACS & NICMOS images of the 
z = 5.83 galaxy, SBM03#1. Right: the Spitzer and HST wavebands straddle the 
age-sensitive Balmer/4000 A break, and reveal an underlying old (~ 400 Myr) pop- 
ulation which dominates the stellar mass, ~ 3 x 10 10 M & (Eyles et al. 2005). 

We have recently discovered evidence of an early epoch of star formation (Eyles 
et al. 2005). Some of the spectroscopically-confirmed i'-drops at z ~ 6 show 
evidence of another spectral break in the near-infrared (comparing Spitzer 
images with IRAC at 3.6 - 4.5 fim with HST NICMOS & ACS images at 
0.8 — 1.6 /xm). This age-sensitive Balmer/4000 A break indicates the IR light 
is dominated by stars older than 100 Myr (and most probably 400 Myr old). 
This pushes the formation epoch for these galaxies to z {orm = 7.5 — 13.5. 
Our recent results also indicate that there were at least some galaxies already 
assembled with stellar masses ~ 3 x 10 10 M o (equivalent to 0.2 M* today) 
within the first billion years. The early formation of such systems may have 
played a key role in reionizing the Universe at z ~ 10. 



References 

[1] Becker R. H. et al, 2001, AJ, 122, 2850 

[2] Beckwith S., Somerville R., Stiavelli M., 2003, STScI Newsletter vol 20 
[3] Blanton M. R., et al, 2003, ApJ, 592, 819 
[4] Bouwens R., et al., 2003, ApJ, 595, 589 

[5] Bouwens R., Broadhurst T., Illingworth C, 2003, ApJ, 593, 640 
[6] Bouwens R., et al, 2004, ApJ, 606, L25 

[7] Bunker A. J., Stanway E. R., Ellis R. S., McMahon R. C, McCarthy 

P. J., 2003, MNRAS, 342, L47 
[8] Bunker A., Stanway E., Ellis R., McMahon R., 2004, MNRAS, 355, 374 
[9] Bunker A. J., Stanway E. R., 2004, |aBtro-ph/0407562| 



10 



[10] Connolly A. J., Szalay A. S., Dickinson M., Subbarao M. U., Brunner 

R. J., 1997, ApJ, 486, Lll 
[11] Coleman G. D., Wu C.-C, Weedman D. W., 1980, ApJS, 43, 393 
[12] Dickinson M., et al., 2004, ApJ, 600, L99 
[13] Dijkestra M., Haiman Z., Loeb A., 2004, ApJ, 613, 646 
[14] Eyles L., Bunker A., Stanway E., Lacy M., Ellis R., Doherty M., 2005, 

MNRAS m press, astro-ph/0502385 
[15] Fan X., et al., 2002, AJ, 123, 1247 

[16] Fernandez-Soto A., Lanzetta K. M., Chen H.-W., MNRAS, 342, 1215 
[17] Fontana A., et al, 2003, ApJ, 587, 544 

[18] Giavalisco M., Dickinson M., 2003, in "The Mass of Galaxies at Low and 

High Redshift" Proceedings of the ESO Workshop, p324 
[19] Giavalisco M., et al, 2004, ApJ, 600, L103 
[20] Gnedin N. Y., Ostriker J. P., 1997, ApJ, 486, 581 
[21] Gunn J. E., Peterson B. A., 1965, ApJ, 142, 1633 
[22] Iwata L, et al, 2003, PAS J, 55, 415 
[23] Kogut A., et al, 2003, ApJS, 148, 161 

[24] Lanzetta K. M., Yahata N., Pascarelle S., Chen H.-W., Fernandez-Soto 

A., 2002, ApJ, 570, 492 
[25] Lilly S. J., Tresse L., Hammer F., Crampton D., Le Fevre O., 1995, ApJ, 

455, 108 

[26] Madau P., Ferguson H. C, Dickinson M. E., Giavalisco M., Steidel C. C, 

Fruchter A., 1996, MNRAS, 283, 1388 
[27] Madau P., Pozzetti L., Dickinson M., 1998, ApJ, 498, 106 
[28] Madau P., Haardt F., Rees M., 1999, ApJ, 514, 648 
[29] Malhotra S., et al., 2005, ApJ, 626, 666 
[30] Pirzkal N., et al, 2005, ApK, 622, 319 
[31] Ricotti M., 2002, MNRAS, 336, L33 
[32] Salpeter E. E., 1955, ApJ, 121, 161 
[33] Spergel, D. N., et al, 2003, ApJS, 148, 175 

[34] Stanway E. R., Bunker A. J., McMahon R. G., 2003, MNRAS, 342, 439 
[35] Stanway E. R., et al, 2004a, ApJ, 604, L13 

[36] Stanway E. R., Bunker A. J., McMahon R. G., Ellis R. S., Treu T., 

McCarthy P. J., 2004b, ApJ, 607, 704 
[37] Stanway E. R., McMahon R. G., Bunker A. J., 2005, MNRAS, 359, 1184 
[38] Stiavelli M., Fall S. M., Panagia N., 2004, ApJ, 610, LI 
[39] Steidel C. C, Pettini M., Hamilton D., 1995, AJ, 110, 2519 
[40] Steidel C. C, Giavalisco M., Pettini M., Dickinson M. E., Adelberger 

K. L., 1996, ApJ, 462, L17 
[41] Steidel C. C, Adelberger K. L., Giavalisco M., Dickinson M. E., Pettini 

M., 1999, ApJ, 519, 1 
[42] Steidel C. C, Pettini M., Adelberger K. L., 2001, ApJ, 546, 665 
[43] Williams R. E., et al., 1996, AJ, 112, 1335 
[44] Yan H., Windhorst R. A., Cohen S., 2003, ApJ, 585, L93 
[45] Yan H., Windhorst R. A., 2004, ApJ, 600, LI 



11 



